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ABSTRACT 

The impact of radiation pressure on the dynamics of the gas in the vicinity 
of young stellar clusters is thoroughly discussed. The radiation over the ther- 
mal/ram pressure ratio time evolution is calculated explicitely and the crucial 
role of the cluster mechanical power and of the strong time evolution of the 
ionizing photon flux and of the bolometric luminosity of the exciting cluster is 
stressed. It is shown that radiation has only a narrow window of opportunity to 
dominate the wind-driven shell dynamics. This may occur only at early stages 
of the bubble evolution and if the shell expands into a dusty and/or a very dense 
proto-cluster medium. The impact of radiation pressure on the wind-driven shell 
becomes always negligible after about 3 Myr. Finally, the wind-driven model 
results allow one to compare the model predictions with the distribution of ther- 
mal pressure derived from X-ray observations. The shape of the thermal pressure 
profile allows then to distinguish between the energy and the momentum domi- 
nated regimes of expansion and thus conclude whether radiative losses of energy 
or the leakage of hot gas from the bubble interior have been significant during 
the bubble evolution. 

Subject headings: galaxies: star clusters: general - HII regions - hydrodynamics 
- ISM: bubbles - ISM: kinematics and dynamics 

1. Introduction 

Several driving mechanisms have been proposed to explain the origin and evolution 
of large expanding structures discovered in deep Halpha images and HI surveys of the 
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Milky Way and other local group galaxies. The most thoroughly discussed model involves 
the thermalization of the kinetic en ergy released by stellar winds and supernovae explo- 
sions inside a compact star cluster (IMcCray &: Kafatod 119871 ; iMac Low &: McCrayl Il988l ). 
While this model is broadly consistent with observations of many bubbles in the visi- 
ble and X- ray emission, several others driving forces, such as radiation pressure from the 
field stars (lElmegreen fc Chiang] Il982l) or the encounter of a high velocity cloud with the 
galactic disk (iTenorio-Tagl e 1981), were also discussed in the literature (see for a review 
Tenorio-Tagle fc Bodenheimerl " 1988 ). The impact of radiation pressure on the dynamics of 
giant HII regions has been discussed recently by several authors who claimed that radiation 
pressure may play a crucial rol e in shaping the in terstellar medium (ISM) around young 
massive clusters. For example, lLopez et al.l ( 1201 ll ) compared the thermal pressure in the 
X-ray emitting plasma with the flux of radiation energy in the 30 Dor reg ion and claimed 
that t he hot gas pressure is generally weak and not dynamically important. lYeh &: Matzner 
( 120121 ) reached the same conclusion by making use the ionization parameter technique and 
claimed that in many individual and averaged over the galactic scale targ ets the shocked 



wind pressure cannot be large compared to radiation pressure. However iPellegrini et al. 



( 120111 ) measured the ionization parameter across the 30 Dor region and concluded that ra- 
diation pressure does not currently play a major role in star forming regions such as 30 Dor, 
although it may have been an important factor during the early evolution. The revision of 
the classical model is thus required for better understanding of such controversial results. 

The development of HII regions begins right from the formation of the stellar cluster, 
when the ample supply of ioniz ing photons leads to a supersonic weak R ionization front 
( lKahnlll954j ; lTenorio-Tagldll976l ) . supersonic with respect to the neutral gas ahead and also 
supersonic with respect to the ionized gas behind it. The ionization front rushes then through 
the surrounding gas leaving it warm (T ~ 10 4 K) and ionized but dynamically unperturbed. 
This situation prevails throughout the formation phase, until the ionization front reaches 
the Stromgreen radius. The ionization front becomes then of D type with a strong shock 
progressing ahead of it into the surrounding gas. From then onwards and through the whole 
of the expansion phase, only a small fraction (less than 1%) of the ionizing photon flux will 
impinge on the ionization front and exert a pressure on the surrounding shell of swept up 
matter. For such a reason, in the classic HII region model one cannot use the full bolometric 
luminosity of the central clu ster to estimate the s t rength of radiative pressure during the 
expansion of an HII region. iKrumholz fc Matzner! (120091 ) assumed that radiation pressure 
causes the ionized gas to pile up into a thin outer shell and then concluded that while 
radiation pressure is generally not important for HII regions around a single or a handful of 
massive stars , it dominates the dynamics of giant HII regions driven by massive star clusters. 
Draind ( 1201 ll ) calculated the density distribution in dusty HII regions with radiation pressure 
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in the isothermal static approximation and showed that radiation pressure indeed produces 
a density gradient. However, the deviation from the homogeneous distribution is moderately 
small unless the parameter A = Qm^mi-, where Q 49 is the number of ionizing photons in 
units of 10 49 s -1 and uhii is the mean density in the HII region, is large (Ao > 10 4 cm -3 ). 

A more realistic and also very different situation arises if one considers also the mechan- 
ical power of the exciting cluster and the strong evolution that the ionizing photon flux, the 
bolometric and the mechanical luminosities suffer after the first supernova explosion. Soon 
after the cluster wind interacts with the ionized surrounding gas, a strong shock wave begins 
to form a secondary wind-driven shell with the surrounding swept up matter. As the situ- 
ation evolves, the pressure in the inner shell becomes larger than that of the free wind and 
thus a second, reverse shock is established in order to equalize and maintain an even pressure 
between the shocked wind and the swept up ionized gas. The wind-driven shell sooner or 
latter cools down and begins to recombine, depending on the value of the background den- 
sity. Such recombinations deplete photons from the outer HII region and cause the ionization 
front to supersonically recede towards the cluster to finally become trapped within the inner 
wind-driven shell. The contribution of the star cluster UV radiation to the expansion of the 
original HII region becomes then insignificant. From then onwards the wind-driven shock 
travels through a warm neutral or molecular gas whose sound speed is much smaller than 
that in the ionized interstellar medium. This per mits the wind-driven bubble to survive for 



a much longer time and lead s to the regime which ICapriotti &: Kozminskil (120011 ) named the 



"best case for wind" model. ICapriotti &: Kozminskil (120011 ) concluded that in this case the 



dynamical evolution of the nebula is dictated by the stellar wind, but noted that the "best 
case for wind" scenario requires the product Lq^uism, where L 9 i 2 is the luminosity of the 
ionizing source in the Lyman continuum and uism is the density of the ambient medium, 
to be large. In the case of a single massive star, this requires a very high density ambient 
medium as the bolometric luminosities are rather small even in the case of the most mas- 



sive O-stars. lArthurl ( 120121 ) took into consideration the effects of mass loading and thermal 
conduction and developed a detailed numerical model for the stellar wind bubble around 
O'Ori C, the main exiting star of the Orion nebula. The calculations confirmed the "best 
case for wind" scenario and demonstrated that radiative heating makes the ionized part of 
the wind-driven shell thick. However, neither the ionizing source time evolution nor the 
radiation pressure effects have been taken into consideration in these calculations. 

Here we consider the impact that radiation pressure provides on bubbles inflated by 
winds driven by massive star clusters. In such a case, the ionizing (L912) and the bolometric 
luminosity of a typical coeval cluster fall down rapidly after the first supernova explosion, 
whereas the mechanical luminosity of the cluster is maintained at a nearly constant level. One 
can notice this in Figure 1 which presents the results from the Starburst99 synthetic model 



-4- 



( ILeitherer et al.lll999f ) for a 1O 6 M coeval stellar cluster with stellar metallicity Z = 0AZ & 
and standard Kroupa initial mass function. We also display on the right-hand panel of this 
figure the L bo[ over L mech and the L 912 over L mech ratios as these quantities are central to 
our further discussion and lead us to conclude that radiation pressure has a poor impact on 
the expansion of the wind-driven shell and thus on the global dynamics of giant HII regions. 




Fig. 1. — The star cluster luminosities time evolution. The bolometric, Lyman continuum 
and mechanical luminosities of a 10 6 M coeval stellar cluster with stellar metallicity Z = 
0AZ & and standard Kroupa initial mass function as predicted by the Starburst99, v6.0.3. 
The solid and dashed lines on the right-hand panel display the ratio of L& j over L mec h and 
L912 over L mec h also plotted function of time. 

The paper is organized as follows: the inner structure, the impact from radiative cooling 
and the time evolution of star cluster driven bubbles in the two possible (energy-dominated 
and momentum-dominated) hydrodynamic regimes are discussed in section 2. The impact 
of radiation pressure on the wind-driven shell is thoroughly discussed in section 3, where we 
calculate the radiation over the dynamical pressure ratio and discuss how this ratio evolves 
with time. Section 4 presents the distribution of thermal pressure obtained from the wind- 
driven model and derives its appearance when projected onto the plane of the sky. We 
also show in this section that the shape of the projected thermal pressure profile allows 
one to distinguish between bubbles evolving in the energy and in the momentum-dominated 
regimes. Our results are compared with observations and other theoretical models in section 
5 and a summary of our major conclusions is given in section 6. 
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The inner structure and evolution of star cluster driven bubbles 



The thermalization of the kinetic energy supplied by stellar winds and supernovae ex- 
plosions within the exciting cluster leads to a large central pressure that causes the exit 



of the reinserted matter as a star cluster wi nd (j Chevalier fc Cleg's 1985 ; Canto et al 



2000 



Tenorio-Tagle et al.l 120071 ; ISilich et al.l I2011T). Several physically distinct regions are then 
formed in the recently ionized gas (j Weaver et al.lll977l ). The innermost (free wind) zone is 
occupied by the high temperature reinserted plasma whose density, temperature and pressure 
asymptotically fall as r~ 2 , r~ 4 / 3 and r~ 10 / 3 while the wind velocity monotonically increases 
to re ach its terminal va l ue, ~ 3 1 / 2 c n , w here c is the sound speed at the star cluster center 



(e.g. 



Canto et al.ll2000i ; ISilich et al.ll201ll ). The encounter of such high velocity outflow with 
the ambient, uniform density ionized medium, leads to the formation of the leading and re- 
verse shocks. The leading shock sweeps up the ionized gas and compresses it into a secondary 
inner dense shell. The reverse shock decelerates and re-heats the free wind matter what re- 
sults in a high thermal pressure in the shocked wind region. The thermal pressure in this zone 
is almost homogeneous as the temperature and the sound speed in the shocked wind plasma 



Chu et al. 


1995; 


Canto et al. 


2000; 


Stevens & Hartwell 



2003 



Silich et al.l 120051 ) due to the large temperature in these zones (T ~ 10 6 K - 10 7 K). 
Thus the interior of the secondary, wind-driven shell is usually transparent to the ionizing 
radiation from the central cluster. 



2.1. Energy dominated wind-driven bubbles 



The theo ry of interste 



developed b y I Weaver et al. 



l ar wind-driven bubbles with a con s tant energy input rate was 



(1977). 



Mac Low fc McCravi (119881 ). iKoo & McKeel (119921 ) (see 



for a review iBisnovatyi-Kogan fc Silichl Il995l ). In this case the radius, Rb, the expansion 
velocity of the wind-driven shell, and the thermal pressure inside the bubble volume, Pb, 
are: 



Rb 



375(7 - l)£mecfe 

_28(9 7 -4)7r P/5M _ 
3 Rb 



1/5 



t 3 / 5 , 



Pb 



7 1/3 
' PlSM 



3(7 — l)L mec h 
28(97 - 4)7^2 



-1 2/3 



(1) 
(2) 
(3) 
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where t is the dynamical time, Pism is the ambient gas density and L mec h is the mechanical 
luminosity of the central cluster. One can find the location of the reverse shock from the 
relation P b = P ram , where P ra m = PwV^ is the free wind ram pressure and the density 
in the wind, /?,„, is calculated at the reverse shock position: p w = L mech /27iRj iS V^ (e.g. 



Weaver et al.lll977f ): 



T \ V2 

lj mech 



The wind-driven shell of swept up interstellar matter is initially adiabatic and hot and 
thus transparent to the io nizing radiation. However, it cools down in a short time scale 



(Mac Low & McCravill988f ): 

n = (2.3 x ltf)Zj$£nj2£l%»yr, (5) 

where ujsm and Zjsm are the atomic number density and metallicity in the surrounding 
medium, L 38 is the star cluster mechanical luminosity in units of 10 38 erg s _1 and 7 = 5/3. 
As soon as the wind-driven shell cools down, it recombines and begins to absorb ionizing 
photons emerging from the central cluster. This reduces the number of Lyman continuum 
photons in the outer HII region. From then onwards, the ionization front detaches from 
the outer shell and moves back towards the wind-driven shell reaching it when this grows 
thick enough as to absorb all ionizing photons from th e central cluste r. One can obtain the 



characteristic trapping time, r trap , from the condition (jComeronl 119971 ): 



N 912 = AnR 2 AR/3n 2 . (6) 

where N 9 i 2 is the number of Lyman continuum photons, = 2.59 x 10~ 13 cm 3 s _1 is the 
recombination coefficient to all but the ground level. We assume that thermal pressure in 
the ionized shell is uniform and equal to that in the shocked wind region (see Arthur, 2012). 
The density of the ionized shell then is n s = Pb/piC 2 HII , where pi = 14mH/^ is the mean 
mass per ion, chii = (kT HII / p t ) l l 2 is the isothermal speed of sound in the ionized gas and 
p t = 14/ '23m h is the mean mass per particle in the completely ionized plasma with one 
helium atom per every ten hydrogen atoms. The shell thickness, AR, is calculated from 
mass conservation: 

n/5M^ (7) 
n s 3 K ' 

Equations CQ), (ED and © then yield: 

(9 7 -4) p,c 2 HII N 9 



-HII 1 " 912 

1 trap 



5(7 - 1) PniSM L m ech 

During the forthcoming evolution the HII region is bounded by the wind-driven shell. The 
leading shock then remains strong even when the shell velocity drops below chii as the 
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leading shock Mach number must be calculated with respect to the sound speed in the 
neutral ambient medium. Note, that the ioniz ing radiation may b e trapped within the shell 



at much earlier times if the shell is dusty (see iMurray et al.l 1201 if ) . 



Similarly, one can find the relative thickness of the ionized shell AR/Rb before and after 
the trapping time r trap . For t < T trap the shell is completely photoionized and the value of 
AR/R b is restricted by the swept up mass M sheU = 47iRlp ISM /3 while for t > r trap the shell 
is partially ionized and the thickness of the ionized layer depends on the number of incident 
ionizing photons A^ 912 : 

AR/R b = ^i^ism t < Ttr ^ (g) 



6 



AR/Rb = %fRT 1 > Ttrap - (10) 



Figure 2 shows the wind-driven shell evolution onto an ambient medium with densities 
u ism — 1 cm -3 and n IS M = 1000 cm -3 (solid and dashed lines, respectively). The selected 
mechanical luminosity of the driving cluster is L mec h = 10 40 erg s -1 in both cases, which is the 
average value for a young 10 6 clust er with a standard Kroupa initial mass function and 
O.4Z metallicity (ILeitherer et al. Il999h whereas the number of ionizing photons changes with 



time as is predicted by the evolutionary synthetic models (see Figure 1). It will be shown 
below that this approximation does not affect our major results significantly because the 
radiative over dynamical pressure ratio depends mainly on the bolometric over mechanical 
luminosity ratio and only weakly on the star cluster mechanical luminosity itself (see section 
3). In the lower density case the wind-driven shell cools after t\ ~ 10 5 yr and grows thick 
enough to absorb all ionizing photons (~ 10 53 s^ 1 ) at about T tra p ~ 3.4 Myr (see Figure 
2). At that time the radius of the shell reaches about 360 pc and its expansion velocity has 
dropped to about 60 km s _1 . Also, the number of Lyman continuum photons emerging from 
the central cluster begins to decrease due to the explosion of the most massive stars (see 
Figure 1). 



2.2. Momentum dominated wind-driven bubbles 

It is likely that young stellar clusters are embedded initially into dense parental clouds. 
Therefore Figure 2 also presents the wind-driven bubble evolution into a dense ambient 
medium with uism = 1000 cm -3 (dashed lines). In this case the shell absorbs all ionizing 
photons as soon as it grows enough to encompass all the cluster. Thus, in the high density 
environment the shell cools down and traps the ionizing radiation from the driving cluster 
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Fig. 2. — The wind-driven bubble evolution. Panels a and b present the radius and the 
velocity of the shell when it expands into an ambient medium with densities: ujsm — 1 cm -3 
and tiism = 1000 cm -3 (solid and dashed lines) respectively. The mechanical luminosity of 
the driving cluster is the same in both cases: L mec h = 10 40 erg s _1 . The metallicity of the 
ambient gas is Z IS m = O.4Z , the rate of ionizzing photons was calculated by making use of 
Starburst99 synthesis model. The arrows mark the characteristic trapping time, T tTapi in the 
lower density case and the characteristic time scale for the shocked wind gas to cool down, 
r 2 , (see sections 2.1 and 2.2 for more details) in the larger density case, respectively. 



very rapidly, much faster than when it expands into a typical ISM with ujsm — 1 cm" 3 . 
Besides this, in the denser environment the shocked wind region may lose a significant fraction 
of the deposited energy in a short time scale due to strong radiati ye cooling in the conductio n 



zone dominated by mass evaporated from the wind-driven shell (IMac Low & McCraylll988l ): 

r 2 = (1.6 x l^)Zj!l{ 2 \f s [) l L% ll yr. (11) 

In the case under consideration this occurs at r 2 ~ 1.6 Myr, when the shell radius reaches 
about 60 pc and the expansion velocity is ~ 20 km s _1 . Without a pressure support, the 
reverse shock moves rapidly towards the shell and the further expansion is then due to the 
direct impact of the free wind momentum. This may also be the case if a collisionless reverse 
shock at Rrs < Rh is not formed (see the discussion of the reverse shock conditions in 



Capriotti fc Kozminskil 120011 ) . One can neglect this transition phase and compute the shell 
dynamics in the momentum dominated regime by writing down the momentum equation, 
assuming that the transition to the momentum dominated regime occurs at t = r 2 : 



M sh(t)^ = M sh (T 2 )V b (r 2 ) + j M^Kodt', (12) 



where M sh {t) = 4nRl{t)p ISM /3 is the mass of the shell at time t and M w = 2L mech /V^ is 
the rate of mass deposition due to the wind. 
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The integration of equation ffT2]) yields: 



Rh — Rr 



nVoaPlSMRi 



12 QL mec f l T2 \ t 
5 irVoopisMRi) t 2 



1/4 



3L 



mech 



ZkVooPjsmRI 

Lmech 



t + 



3 i£ 3 L mec j l T2 



5 r 2 2 tiVooPism J 



(13) 
(14) 
(15) 



where i? c is the radius of the shell at the time when the transition to the momentum domi- 
nated regime occurs: R c = Rb{j 2 ) and Rb{j 2 ) is calculated by means of equation (pQ). Note, 
that if i? c — ?■ and r 2 — > 0, equation (fT3|) is reduced to a well known relation R b ~ t 1 / 2 (e.g. 
equation 3.1 in lKoo k McKeelll992h . 



Note that in the momentum-dominated regime the value of the thermal pressure in the 
wind-driven shell is equal to the wind ram pressure at the shell inner edge. The relative 
thickness of the ionized layer then is: 



AR/R b 



KP-V^HjjRbNgu 



(3L 2 



(16) 



mech 



If thermal cond uction and mass evaporation from the cold shell are inhibited by magnetic 
fields (jSpitzerl Il962l ) . the radiative losses of energy from the shocked wind region are much 
smaller. In this case the ion number density and temperature in the shocked wind region 
are: 



^PiV^Rl ' 



3Mt 
ArcpiRl 

Pt Pb 

Pi kn b 2(97 - 4) k 



5(7 - l)PtV c 



2 

oo 



;i7) 



;is) 



Thus temperature in the shocked wind region in this case depends only on the wind terminal 
speed and does not change with time. If Voo = 1000 km s -1 , TJ, w 1.1 x 10 7 K. 

Radiative cooling begins to reduce the thermal pressure in the shocked wind zone and 
affect the bubble dynamics only when radiative losses of energy from the bubble interior 
Q = 47™ 2 A(T fe , Z b )Rl/3 exceed the energy input rate from the central cluster: Q > L mech . 
This occurs if the evolutionary time t grows larger than 



t 3 



4 

oc 



V3A(T b ,Z 6 



7T 



J mech 



375(7 - 1) 
28(97 - 4)p 



ISM 



(19) 
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where Zb is the gas metallicity in the shocked wind region and A(Tb, Zb) is the cooling 
function. One can obtain from equation f fl9|) that r 3 is much longer than the characteristic 
life-time of giant HII regions (~ 10 Myr). This implies that the radiative losses of energy from 
the shocked wind region in this case are insignificant and the wind-driven bubble expands 
in the energy dominated regime during the whole evolution of the HII region. 

The reverse shock may also be closer to the wind-driven shell than one would expect from 
equation (TJJ, if the hot, shocked win d plasma escapes fr om the bubble into the surrounding 



low density medium, as suggested in lLopez et al.l (1201 lh . In this case the expansion velocity 



of the wind-driven shell differs from that predicted by the energy dominated model (equations 
CD -[31) and may be more similar to what the momentum dominated model predicts. 



3. The impact of radiation pressure on the wind-driven shell 

The radiation pressure on the wind-driven shell is: 

Prad = ftrapLbol/^TlcRl, (20) 

where f trap is the fraction of the bolometric luminosity absorbed within the shell, c is the 
speed of light and Rb is the radius of the shell. In the energy dominated regime equations 
(U and (T20D yield: 

Prad ftrap ( Lbol \ ( Voo \ / Rrs 



■ ( 21 ) 

A 2 \L 

mech J \ C« , I V Rb 

One can note from relation (pTj) that the small over c ratio and the fact that Rrs is 
usually much smaller than Rb, decrease significantly e from what one would expect during 
the early bubble evolution when the bolometric luminosity of the driving cluster exceeds the 
mechanical energy input rate by a factor of ~ 1000 (see Figure 1). Hereafter we assume in 
our calculations that = 1000 km s™ 1 . 

The P ra d over Pb ratio may be presented in a different form which shows explicitely how 
it evolves with time, if one makes use of equation Q instead of equation (TJJ: 



rad 



ftrap 
287TC 



28(97 - 4)vr 
3(7 - 1) 



4/5 / j \ 1/5 



Lmecn j ^ooi {22) 

25piSM/ L 

mech 



Equation ( 1221) shows that the ratio of the bolometric over mechanical luminosity is the 
major factor which decides which of the two driving forces, radiation or thermal pressure, 
dominates the shell dynamics and that the P ra d over Pb ratio has only a weak dependence 
on the ambient gas density and on the star cluster mechanical luminosity. 
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The relation between P rac i and Pb will be different if the shell expands in the momentum 
dominated regime. In this case, one has to compare the radiation and the free wind ram 
pressure at the inner edge of the shell: 



Prad ftrap ^oo Lbol 



Pram 2 C i-'mech 



(23) 



Note, that relationship ( 12 3 p is universal as it does not depend on the ambient gas density 
and the Lboi over L mec h ratio is a universal function of time for all clusters with a given initial 
mass function and metallicity. 

Figure 3 presents the P ra d over Pb ratio time evolution for the two cases here considered. 
As mentioned above, it was assumed that the L bo i over L mech changes with time as predicted 
by the evolutionary synthesis model Starburst99 but a constant value of the mechanical 
luminosity, L mec h = 10 40 erg s -1 , was used in the calculations as e depends on the value of 
Lmech itself only in the energy dominated regime and even then this dependence is weak (see 
equations [22] and [23]). Numerical solutions with Starburst99-generated input mechanical 



luminosity were used by iDopita et al.l (120051 ) in their discussion of the synthetic spectral 



energy distribution in starbursts. 

Note that Pb in Figure 3 is the thermal pressure in the shocked wind region in the energy 
dominated case and Pb = P ra m m the momentum dominated regime, respectively. In the 
low density case (solid line) the P Ta d over Pb ratio is larger than unity only at very early 
times it < 2 x 10 5 yr) and only if the shell is dusty and optically thick to the bolometric 
luminosity of the central cluster {ftrap = 1)- Otherwise it should grow thick enough to absorb 
all ionizing photons only after T tra p = 3.4 Myr, when the contribution of radiation pressure 
to the shell dynamics is already negligible. The P ra d over Pb ratio is small during the initial 
shell evolution if the driving cluster is embedded into a high density medium (see dashed line 
in Figure 3 and equation [22]) . It then suddenly grows to a value of 1.4 (see equation [23]) as 
the gas in the bubble interior cools down, the reverse shock reaches the shell and the further 
bubble evolution proceeds in the momentum dominated regime. However, even in this case 
the P ra d over Pb ratio becomes smaller than one after ~ 2 Myrs and thus radiation has only 
a narrow window of opportunity to dominate the shell dynamics. Note also that the true 
value of ftrap is set by the SED of the incident spectrum and by the optical depth of the shell 
and may be smaller than unity. In this respect Figure 3 presents an upper limit for the P ra d 
over Pb ratio. For example, if the shell is thick to the Lyman continuum only, the P ra d over 
Pb ratio would be smaller than what is shown in Figure 3 within a factor of 2 - 4, depending 
on the star cluster age (compare the bolometric and the Lyman continuum luminosities in 
Figure 1). On the other hand, the adiabatic wind terminal speed also changes with tim e 



and may be different within a factor of 3 from the selected value (see lWunsch et al.ll201ll ) 




Fig. 3. — The contribution of radiation pressure to the star cluster wind-driven shell dynam- 
ics. Vertical arrows mark the time when the transition to the momentum dominated regime 
occurs (Y2) and the trapping time for the ionizing radiation (r trap ) in the high and low density 
models, respectively. Note that in the low density case (solid line) the wind-driven shell does 
not even trap all ionizing radiation during the early evolution (r trap 3.4 Myr) if it is not 
dusty. The dashed line shows the P ra d over P b ratio in the case when n IS M = 1000 cm -3 . 
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4. Thermal pressure as an indicator of the hydrodynamic regime 

Here we calculate the distribution of thermal pressure inside a wind-driven bubble vol- 
ume and compare it to its proj ection onto the pla ne of the sky as it is done in high resolution 



studies of X-ray emission (e.g. lLopez et al.ll201ll ) 



In order to calculate the distributions of physical quantities (density, temperature, ther- 
mal p ressure, velocity) in the free wind region and localize the reverse shock position, we make 



use of ISilich et al.l (120111 ) star cluster wind driven model, which assumes that stars are expo- 
nentially distributed within the star cluster volume. The model input parameters are: the 
starburst mechanical luminosity, L mech = 10 40 erg s -1 , the characteristic scale length of the 
stellar density distribution, R core = 1 pc, the adiabatic wind terminal speed, Vaoo = 1000 km 
s^ 1 , and the metallicity in the free wind region, which we set to be Zx = 0AZ & . The position 
of the reverse shock at different times was calculated by means of equation which relates 
the thermal pressure in the shocked wind zone and the ram pressure in the free wind region. 
We then compute the density and the temperature in the shocked wind region by means of 
equat ions (fT7|) and (TT8D and calculate the distribution of thermal pressure in the same way, 



as m 



Lopez et al.l (120111 ). Specifically, we calculate the emission measure and the waited 
temperature of the hot X-ray plasma along lines of sight with different impact parameters 
X by integrating the model predicted density and temperature distributions. The emission 
measure and the waited temperature then are: 



-'max 



EM(X) = 2 / n 2 (r)dl, (24) 



T(X) =2EM~ 1 (X) / T(r)n 2 (r)dl, (25) 
Jo 

where I and L max = (R 2 — X 2 ) 1 / 2 are the distance and the path length along the line of 
sight, n{r) and T(r) are calculated at r = (X 2 + I 2 ) 1 / 2 . One can obtain then the plasma 
density and the thermal pressure along lines of sight with an impact parameter X from the 
relations: 

n{X) = [EM{X)/2L max \ 1 ' 2 , P{X) = ^kn{X)T{X)/ix t . (26) 

Note, that the procedure takes into account that lines of sight cross the free wind region if 
the impact parameter X is smaller than the reverse shock radius Rrs- We also take care 
to omit segments along lines of sight where the temperature drops below the X-ray cut-off 
temperature, T cut = 5 x 10 5 K (see Figure 4). 

In the momentum-dominated regime the shocked wind region collapses and the star 
cluster wind impacts directly on the shell. Strong radiative cooling reduces then the temper- 
ature immediately behind the reverse shock leaving the shocked wind matter photoionized 
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and at 10 4 K. Thus, in the momentum-dominated regime the contribution of the shocked 
wind gas to the X-ray emission is negligible. 

The distributions of thermal pressure obtained from the wind-driven bubble model and 
those, obtained by integration along different lines of sight are compared in Figure 4 (upper 
and bottom panels, respectively). The middle panels present the temperature distribu- 
tions, which one can obtain from the wind-driven bubble model. Panels on the left-hand 
and the right-hand sides compare profiles obtained in the low (tiism — 1 cm~ 3 ) and high 
{nisM — 1000 cm~ 3 ) density cases. Solid and dashed lines show the distributions of ther- 
mal pressure and temperature at different evolutionary times: t = 2 Myr and t — 5 Myr, 
respectively. Figure 4 shows that the physical properties of the X-ray plasma derived from 
observations must be taken with care as the distribution of thermal pressure may be sig- 
nificantly distorted and the true central pressure significantly underestimated due to the 
inhomogeneous distribution of plasma inside the wind-driven bubble volume. For example, 
in the energy dominated regime the value of central pressure obtained by integration along 
lines of sight is within a factor of ten smaller than the actual value of the thermal pressure 
at the star cluster center. The integrated profile is significantly shallower and does not show 
a deep gap, which one can notice in the model predicted distribution of thermal pressure. 
The minimum in the integrated distribution of thermal pressure is shifted from the model 
predicted position towards the center and therefore it does not mark the position of the 
reverse shock (compare the top and the bottom panels in the left-hand column of Figure 4). 

One can also note in Figure 4 that the shape of the integrated profiles is very different 
in the energy (low density) and momentum (high density) dominated cases. In the energy 
dominated regime the distribution of thermal pressure has a strong maximum in the center 
and a uniform plateau at larger radii. In the momentum dominated case the integrated profile 
is more similar to what the free wind model predicts (compare the top and the bottom panels 
in the right-hand side of Figure 4). In this case the major difference between the model and 
the integrated profiles occurs in the center, where the model predicted pressure is about ten 
times larger than that obtained by integration, as it also occurs in the energy-dominated 
case. 

Thus, the shape of the thermal pressure profile allows one to distinguish between the 
two possible hydrodynamic regimes. The thermal pressure profile with a narrow central spike 
and a uniform plateau implies that the wind-driven bubble evolves in the energy dominated 
regime and that the reverse shock stands within the explored volume. A thermal pressure 
that drops continuously with distance from the star cluster center indicates that the shocked 
wind zone has collapsed due to strong radiative cooling and the wind-driven bubble evolves 
in the momentum dominated regime. The reverse shock accelerates then out from its energy- 
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Fig. 4. — The model predicted vs integrated thermal pressure profiles. The upper and 
middle panels display the model predicted distributions of thermal pressure and temperature, 
respectively. The thermal pressure profiles obtained by integration of the model predicted 
quantities along different lines of sight (which one might expect to obtain from the observed 
X-ray emission) are shown in the bottom panels. The left-hand and right-hand side columns 
present the results for the low and high density (iiism — 1 cm -3 and tijsm — 1000 cm -3 ) 
models, respectively. The solid and dashed lines display profiles at different evolutionary 
times: t — 2 Myr and t — 5 Myr. The horizontal solid lines in the middle panels display the 
X-ray cut-off temperature. 
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dominated position and finally the star cluster wind ends up impacting directly on the outer 
shell. 



5. Comparison to other results 



The impact f rom radiation pressure on th e dynamics of giant HII regions has recently 
been discussed by iKrumholz &: Matznerl (120091 ) who claimed that radiation pressure might 
dominate the dynamics of HII regions around massive star clusters, whereas we found that 
radiation has only a narrow window of opportunity to dominate the giant HII regions global 
dynamics. 



The key difference between our results is that IKrumholz fc Matznerl (120091 ) did not take 
into consideration the mechanical power of the exciting cluster and the time evolution of 
the star cluster bolometric luminosity. Thus, the radiation force in their and our models 
is compared with different pressures: with thermal pressure in the ionized gas behind the 
ionization front in the Krumholz & Matzner model and with thermal/ram pressure in the 
shocked/free wind zone in the wind-driven bubble model. Another significant difference 
is that we consider the time evolution of the exciting cluster explicitely that restricts the 
radiative feedback as the bolometric luminosity and the flux of ionizing radiation drop rapidly 
in massive coeval clusters with a normal IMF and ages larger than 3 Myrs. 



The i mpact of r adiati ve heating on the bubble dynamics has been discussed by lCapriotti fc Kozminski 



(120011 ) and lArthurl (120121) who have considered b ubbles driven by a single massive star. As 
it was stressed by ICapriotti fc Kozminskil ( 120011 ). the real distinction between the ionizing 
radiation and the star cluster wind feedback to the ambient medium is that the character- 
istic lifetime of the star is much longer than the wind-blown bubble stalling time, the time 
when the bubble expansion velocity drops to the sound speed value in the ambient medium 
and the wind-induced shock dissipates. This is, however, not true in the case when bub- 
bles are driven by massive star clusters as the characteristic stalling time increases with the 
mechanical luminosity of the driving cluster: 

1/2 

(27) 



TstalLED ~ 2.7c 



-5/2 
ISM 



(7 - ^)L me ch 



_7(97 - 4:)np IS M_ 



T stalLMD 



0.25Cjg M 



3L 



mech 



nl/2 



where T s t a ii,ED and t stall, md are stalling times in the energy and momentum-dominated 
regimes, respectively, cism is the sound speed in the ambient medium and we have used 
an asymptotic >> R c ) relation for a momentum-dominated bubble expansion velocity 
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(see equation [Ml) . In the case of massive star clusters the characteristic trapping time Tt rav 
is short compared with both, the wind-driven bubble stalling time, r sta u, and the charac- 
teristic lifetime of the HII region, r HII w 10 Myr (e.g. in our models T sta u,MD ~ 13 Myr 
and T st au,ED > 40 Myr if T IS m = 1000K). This implies that in this case wind-blown bubbles 
expand for a rather long time in the "best case for wind" regime, when the wind-driven shell 
traps all ionizing photons from the central source and the leading shock Mach number must 
be calculated with respect to a low sound speed in the neutral ambient medium. 



In the calculations by I Arthur! (120121 ) the ionized shell grows thick due to radiative heating 
what affects the wind-driven bubble dynamics because it reduces the volume occupied by 
the shocked wind gas and thus enhances the thermal pressure in the shocked wind zone and 
behind the leading shock front. This effect, however, strongly depends on the mechanical 
luminosity of the driving cluster because the thickness of the ionized s hell is smaller in systems 
with a higher internal pressure (compare Figures 2-3 and 4-5 in lArthurl 120121 ). Figure 5 
shows the relative thickness of the photoionized shell calculated by means of equations (Q - 
(flQj) and (Ti~6l) . As one can see, in the case of a 1O 6 M coeval cluster, the thickness of the 
ionized shell does not exceed ~ 10% of the wind-blown bubble radius what implies that in 
this case radiative heating does not affect significantly the bubble dynamics. Equations (J2J 
- ffTOl and f|T6|) predict that the AR/R b ratio grows larger for low mass/energetic clusters. 
Thus, it is expected that the HII gas filling factor, fan, has to be larger in the HII regions 
developed by the less massive c lusters as it is the case in the Carina nebula, where fun ~ 1 
(IHarper-Clark &: Murray! 120091 ) and in the 30 Dor regio n, where the filling factor of the 
ionized gas is very small: fun ~ 0.03 (jMills et al.lll978l ). The effects of radiative heating 
must be also reduced if the wind-driven shells are dusty as the dust competes with the gas 
for ionizing photons what reduces the number of the incident photons heating the shell. 



Harper-Clark fc Murray! (120091 ) discussed the dynamics and the X-ray emission from the 



Carina Nebula bubble and claimed that the leakage of the shocked wind gas through the holes 
in the expanding s hell may solve a long t ime standing problem of the shell growth descipancy 
what implies that IWeaver et all (119771 ) model oft en overpredicts the sizes and expansion 
velocities of the wind-driven shells (see IOeylll996l ). The leakage gas model, however, does 
not work in this case as the reverse shock radius calculated by means of equation (jlj) with 
the v alues of thermal pressure , Tr 16 cluster mechanical luminosity and mass-loss rate taken 
from lHarper- Clark fc Murray! (120091 ): P b = 2 x 10~ 10 dynes cm -2 , L mech = 3.5 x 10 38 erg 
s _1 and M to t = 1.1 x 10 -3 M yr _1 is comparable or even exceeds the observed size of 
the Carina bubble (Rrs ~ 17 pc while the observed values of the Carina bubble radius 
range from 10 pc to 20 pc). The leakage bubble model then cannot be used in order to 
evaluate the X-ray luminosity and brightness profile from the Carina Nebula as it does 
not account for the free wind zone (it is assumed that Rrs << Rb) and thus neglects the 
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Fig. 5. — The relative thickness of the ionized shell. Solid and dashed lines present the 
AR/R b ratio evolution in the case when the driving cluster is embedded into a low (1 cm" 3 ) 
and high (1000 cm" 3 ) density ambient medium, respectively. 
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contribution of the free wind region to the X -ray emission. Note also that the apparent 
discrepancy between the IWeaver et al.l (119771 ) model predictions and the observed bubble 
radii and expansion velocities may b e understood if bubbl es evolve from flattened parental 
clouds and are observed face-on (see 



Silich fc Franco! Il999h 



The geometrical ionization parameter U = Ngi2/4ircr 2 nHii = N^k^iTi/ ' ^nxc^tPbi where 
c is the speed of light, reaches the maximum value in the momentum-dominated, high density 
model: loq(U m n . T ) ~ This value is smaller than the max i mum value obtained by 

Yeh fc Matznerl (120121 ). log{U max ) ~ —1, who used the iDraind (120111 ) static model and 



assumed that the density and thus thermal pressure are zero at the inner edge of the strongly 
illuminated HII region (see A ppendix B in the paper). This may be the case in a static, 
pressure confined HII regions (lDrainell201ll ). but is not valid in the case of the wind-driven 
shell where the thermal pressure at the inner edge of the shell must be equal to that in the 
shocked wind zone even if the density distribution in the ionized shell is not homogeneous 
due to strong radiation pressure. 

A flat plateau and a steep central spike predicted by the energy-dominated wind-driven 
bubble model agree surprisingly well w ith the shape of the thermal pressure profile of 30 Dor 
region obtained by [Lopez et al.l (120111 ). This allows one to rule out the momentum domi- 
nated regime and conclude that the leakage of the hot gas from the 30 Dor shell is not a 
significant factor as suggested in lLopez et al.l (120111 ) . This conclu sion also agrees w i th the 
deta iled analysis of t he X- ray emission from the 30 Dor region by IWang fc Helfandl (Il99ll ) 
and iTownsley et al.l (120061 ) who claimed that the diffuse X-ray emission is enveloped and 
probably confined by the cooler gas which outlines the classic picture of 30 Dor and extends 
for about 300 pc from north to south. 

Radiation pressure is approximately equal to the gas pressure in the highly ionized ridge 
in th e center of 30 Dor reg ion. Outside of this central zone the P r ad/Pgas ratio drops below 
1/3 (Pellegrini et al.ll201lf) and the ionization parameter U does not show much gradient 
(llndebetouw et al.l 120091 ) . The P ra d/Pgas < 1/3 value is roughly consistent with the wind- 
driven model predictions (see Figure 3), if one bears in mind the complicated history of star 
formation in the 30 Dor starburst with three episodes of star formation o f increasing intensit y 
which took place approximately 5 Myr, 2.5 Myr and about 1.5 Myr ago (jSelman et al.lll999l ). 
The radius of the central arc (~ 10 pc) is comparable to the size of NGC 2070, the major 
driving cluster in the 30 Dor starburst, what requires a more re alistic than a spheric al shell 
illuminated by a central source model as was also stressed by ISnijders et al.l (120071 ) when 
discussing the large values of the ionization parameter {log{U) > —1.53) detected in the HII 
regions around two young (ages less than 3 Myr) stellar clusters in the Antenna. 



Our calculations thus favor the results by iPellegrini et al.l (120111 ) who obtained a much 
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lower radiation over thermal pressure ratio in the 30 Dor region than lLopez et al.l ( 1201 ll ). 
Note also that Pellegrini et al. results l ook in better agree ment with the spectroscopic virial 
mass of the 30 Dor cluster obtained by lBosch et al.l (120091 ) (M vir k9x 10 5 M Q ) which leads, 
in the case of a coeval cluster with a Kroupa IMF and 0AZ Q metallicity, to the maximum 
value of bolomet ric luminosity 4.4 x 10 42 erg s _1 what is about two times smaller than the 
Lb i obtained by lLopez et al.l (120 111 ). 



6. Conclusions 

Here we have explored the contribution of radiative and dynamical (thermal or ram) 
pressure to the dynamics of giant HII regions and found only a narrow window of opportunity 
for radiation pressure to be a dominant factor. We stress the importance of the mechani- 
cal feedback from the exciting cluster and the radiation power time evolution, and derive 
analytic relations which show how the radiation over dynamical pressure ratio evolves with 
time. Careful analysis of the two extreme models which cover conditions ranging from the 
interstellar medium in normal galaxies to those found in dense giant molecular clouds, led 
us to conclude that the dynamical pressure dominates always after about 3 Myr, if one uses 
the standard bubble and star cluster synthesis models. By that time the major factor which 
defines the relative contribution from the two driving forces, the bolometric over the mechan- 
ical luminosity ratio, drops within a factor of ten (see Figure 1) and the impact of radiation 
pressure becomes soon negligible even in the most preferable case. We want to stress that 
the reduction of thermal pressure in the shocked wind zone due to leakage of the hot shocked 
plasma out of a porous shell cannot change this c o nclus ion as may only lead to the displace- 



ment of the reverse shock from the IWeaver et al.l (119771 ) model predicted position, but does 



not affect the distribution of the ram pressure in the free wind zone. In the extreme case this 
may lead to the momentum dominated expansion as it is the case if the shocked wind zone 
collapses due to strong radiative cooling or if the re verse shock is not formed due to the large 



proton mean free path in the shocked wind region (ICapriotti fa Kozminskill200ll ). This leads 
to the maximum value of the P ra d over Ph ratio which is shown in Figure 3 and occurs when 
the wind-driven bubble evolves in the momentum-dominated regime. However, even in this 
"best case for radiation" regime the contribution from radiation pressure becomes small after 
about 3 Myr (see dashed line in Figure 3). We thus conclude that radiation pressure, despite 
being significant during the earlier evolution, in general has a poor impact on the expansion 
of giant shells powered by massive star clusters unless the Starburst99 synthesis model over- 
estimates the star cluster mechanical luminosity significantly or a significant fraction of the 
mechanical energy is lost inside the star cluster volume where the stellar mechanical energy 
is thermalized in random collisions of nearby stellar winds and supernovae ejecta. 
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